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Abstract 


The m asses of excited heavy mesons are studied with sum rules in the 


heavy quark e ective theory. A set of nterpolating currents creating (anni 


hilating) excited heavy m esons w ith arbitrary spin and parity are proposed 


and their properties are discussed. Num erical results at the leading order of 


the O (= g ) expansion are obtained for the lowest doubkts (0* ;1* ) and 
(ar ;27 ). 


Typeset using REV TrX 


P roperties of excited heavy m esons com posed ofa heavy quark and a light antrquark 


are of interest for several reasons. Some of these excited states have been observed in 


experm ents. They w illbe ob £cts of further study in future experm ents w ith B -factories. 


Tn particular, they are useful for tagging in CP experm ents H. Theoretically, the relh- 


tive sm plicity ofthe dynam ical condition m akes t of interest for exploring the intemal 


dynam ics of system s containing a light quark. 


The 1=m 9 expansions are useful for such systems m. isthe heavy quark mass). The 


heavy quark e ective theory HQET) al is a powerful tool. However, | 


for obtaining some 


detailed predictions one needs to combine it with some non-perturbative methods. The 
spectra and decay w idths ofheavy m esn excited states have been studied w ith the 1=m 9 
expansion in the relativistic Bethe Salpeter equations in B,4]. They can also be studied 
with the OCD sum rules n HOET which have been used for ground states of heavy 


m esonsB 6]. 


Asa rst step, in the present article we halus the QCD sum wks in HOET to 


obtain the masses of lowest heavy meson doublets (0* ;1*) and (1* ;2*). 


have been studied with the sum niks in the fullOCD in the literature [ff] 


in the follow ing that using the HQET has some advantage for such probkm s. 


= 


In [8] twas shown that ifthe interaction kemel contains only 


These states 
. We hall se 


and I on the heavy 


quark lne, solutions of the BetheSalbeter equation for heavy mesons of arbitrary pin 


and de nite parity consist of two series of degenerate goin doublets in the in nity quark 


massImitmg ! 1 .Thewave functions pi; where j,P are pin and party and i= 172 


for these two series of states, have the follow ng m ost general om 
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In the above form ulas, Mg is the meson mass, x is the coordinate of the heavy quark, u 
is the coordinate of the light antiquark relative to the heavy quark, v is the velocity of 


is the transverse, 


the heavy meson, u = V u,v = u UV, + = &v, , 4 


sym m etric and traceless polarization tensor. Finally, ‚and ; are Invariant functions of 
u? and uj. 

The assum ptions used to obtain these results can be slightly generalized com pared to 
that explicitly written down in [¥]. W e only need to assume that in them, ! 1 mt 
the kemel is of the form 


G P;p;p) =I Us ©rp;p) + Uy P;p;p)) (3) 


where P is the totalm om entum , p and p° are the relative m om enta of the quarks in the 
iInitialand nal state respectively, and are arbitrary functions of momenta and 

matrices for the light quark which transform as a scalar and a vector respectively, Us 
and Uy are arbitrary invariant functions of momenta. Instead of the free propagator, 


the general form can be used for quark propagators appearing in the B-S equation. The 
calculations leading to above results only need to be slightly changed for this generalcase. 
The restriction to I and in (9) forthe heavy quark is equivaknt to that the interaction 
is independent of the heavy quark pin in the lmitmg ! 1 . This is satis ed by QCD f 


Besides this, the interaction kemel is m ost general. The B-S wave function is not gauge 
invariant so that functions ; and ; are gauge dependent. However, the general form of 
(1) and A) is the sam e for any gauge. Physically, each doublet has a de nite value forthe 
total goin of the light com ponent. Except for the states 0 and 0°, foreach state in the 
series (1) there isa state w th the same } P in the series D) . They correspond to di erent 


1 il 
values of } = j+ > and 7 rå h the follow ing we shall sw itch between the indices i 
and 4. This pattem is in agreem ent w ith the heavy quark symmetry. The above results 


are m odel independent consequences of the heavy quark symmetry. The B-S equation, 


Im QCD theek entary interaction vertex of the quark with the guon is , the heavy quark 


1+ E., ae ite 1+ & 1+ & i 
5 mthemg ! 1 a ne 5 =v ‚the contribution 


2 
ofany Feynm an diagram to the B-S kemel can be reduced to the fom p after integrating out 


1 
propagator is —— 
v k 


the intemalm om entum ofthe diagram . 


being Lorentz covariant, only serves as a convenient tool for nding the om ofthe wave 
function of the state with de nite values for the angularm om enta j, Jı and the parity 
P. For related results, see [9f1d]. The authors of these latter works used directly the 
transform ation properties of the wave function of a com posite systen with the angular 
momenta j and %. Ed] contains results sm ilar to ours, though the author did not deal 
exactly with the B-S wave function. This again m plies that our results on the fom s of 
the B-S wave functions in them, ! 1 Im t are m odel independent. 


The relations beween B-S wave functions and the m atrix elements of the currents 
m ake the above results useful for constructing the appropriate interpolating currents and 
deriving their properties. The wave function i in (1) and (J) can be written in the 
form 


jz i5 el ee a) fı (ug ;u) + Ee wu) : (4) 
W e propose to use 
Jipa 7 &) = hy ix) ER 71 SG) ; (5a) 
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as the interpolating current which creates the corresponding heavy meson, where 


, 1+ & 
h, &) = SoS (x) isthe heavy quark ed in HQET andQ (x) isthatin MUOCD, 


@ 
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where = ° Y ° annihilate the samemeson. The cureents de ned in (5) and (6) have 


nice properties. Let J;P ;ii be the heavy meson state w th the quantum numbers j,P,i 


nthem, ! 1 Imt.Wehave 


OY; 70) H;P i = fp a ane (7) 
where fø; = fp; has the same value for the two states in the same doubt. The formula 
M) can be veri ed as the following. W ithout loss of generality we can choose a gauge 
satisfying A (x)= 0. In such a gauge, it follows from the de nition of the B-S wave 
function jpa («;u) = MI © &)qk+ u)) P;ii that fora current in (6) 
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where =m, M..Usngtheßmsof i (1) and 2) and the relations of the type 
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where S denotes sym m etrization and subtracting the trace tem s in the sets ( ı j) 
and (1 39), []) can be obtained from (d). Only the tem fı contrbutes to fp}. If 
instead of J we use J’ in (6X), the calculation is am ost the same, except that now only 
the f? tem contributes to f;,. 


U sing the equation ofm otion forthe light quark ið q= 0 and that forthe heavy quark 
in HQET iv Dh= 0, one can obtain the relation 


fpi= GE ap ifjp i (10) 


between the coupling constants for the current w ith the extra factor ið and that w thout 
qt. 

Sm ilarly, calculated with the quark and gluon eksnthem, ! 1 Imi the one 
lation functions of J;p ; and Typ 0,0 satisfy 


Yu 30 
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dt & vt) Pu); (11) 


where p; isthe same function for the two states in the sam e doubkt. For verifying [L1]) 


we use the Schw inger gauge 


xA = 0: (12) 


The heavy quark propagator n HOET is 


ihoT h,&)h,(0) Pi= dte vt): (13) 


The kading order Lagrangian in HQET isLo= ihv Dh. Therefore in gaugeld) there 
is no gluon Ine tem inating on the heavy quark Ine. W e thus obtain 
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Z x 
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From Lorentz and translation invariance for gauge invariant quantities, S (y;x) has the 


general fom 


1+ @ P a 7 il+ © 
= Si(Y x)fjyr X) t+ i@aS2:(Y ery. X : 


S (y;x) 


Substituting {17) into (19) brings t into a fmm alm ost identical to the rhss: of (9). The 
form ula [U]) then follows after sim ihr calculations. 


(1) and FL) m ply that two curents in the sets [5q) and BH) with not identicalvalues 


ofj,P,inevermixnthem, ! 1 Im either in the heavy meson or in the uark-gluon 


level and heavy quark symmetry is explicit with these currents. This veri es that they 


are the appropriate interpolating currents for heavy meson states with de nite j,P and 
the light quark momentum „which ae conswednoCD nm. ! 1 Imt. 


The properties for currents m entioned above are mportant for applications to OCD 
sum ries for excited heavy mesons. If one use other currents or study the sum rules 


in ulOCD,ther are in general contributions from two nearby poles corresponding to 
states ofthe same j, P . Their contributionsm ay not be separated correctly. Furthem ore, 
the m ixing of such two states can be calculated within our fom alim by introducing the 
O (=ç ) tem s in the Lagrangian in HOQET. 


For obtainng OCD sum rules in the HOET we study the conelator 
Z 


oa RE dx DI 2... 8; 0) Di: (18) 
From [9), it has the follow ing form at the kading order 
wa) OWS DSA? a'g gal); (19) 
where ! = 2v kand all) Sal) is Independent ofj. 


We shall con ne us to the leading order tem s in this article. The contribution to 
p; (| ) of the pol tem on the hadron side is 


2 
fg id 


2 ! 


On the quark-gluon side of the sum wk, »,,(!) is calculated with the kading order 
Lagrangian Lo in the HQET. The Feynm an diagram s contributing to (3) are shown in 
Fig 1. For the sake of sm plicity we consider the doublets (0* ;1*) and (1*;2"!) in the 
follow ing. 


A ccording to [5d) and BH), there are two possible choices for currents for the doubkt 
(0* ;1" ), either 


Joie 2 = Es hg; (20) 


Jina = BS he er. (21) 


or 


0, + 
Jon 27 Bs hy ( i) Dq; (22) 
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Ina = Zu e3) 


Sm ilarly, there are two possibk choices for the currents for the doubkt (1* ;2*). One is 


Jr 24,2, Zu ; (24) 
1j+ ;l 4 Vv u 3 t t I, 

Ss — 
Y 12 — IE ( i ( 1D, 2 + 2D, 1 2g.1 26 ) . (25) 
Zt v9 et t rt 4 tq: 


A nother choice is obtained by adding a factor iB; to 24) and (29). Note that, without 
the last term in the bracket in (24) the current would coupk also to the 1* state in the 
doublet (0* ;1* ) even mthelmitofn nitemg. 


U sually, the currents w ith the kast number of derivatives are used in the QCD sum 
rule approach. The sum mulesw ih them have better convergence in the high energy region 
and often have better stabilty. However, there is a m otivation for using the currents gJ), 
(24) for the doublet (0° ;1*). In the non-telativistic quark m ode], which usually gives 
correct ordering of energy levels of hadron states, the doublets (0* ;1*) and (1* ;2*) are 
orbital p-wave states which correspond to one derivative in the space wave functions. In 
fact, asm entioned above, the coupling constant for the currents PA), EI) is proportional 
to 2 in (2) and that for the cuments PA), RJ) is proportionalto ı in (1) which are 
the large com ponents in the non-elativistic approxm ation. Therefore, we shall consider 
both the currents Pd), RI) and RJ, RJ) forthe doubt (0* ;1"). 

F irst we consider the doublet (0° ;1"). For the correlator 


Z 
op2= i Are Traun) Pi; (26) 


where dr ; is the current con ugate to 23) ‚thesum rules after the B oreltransfom ation 
is found to be 
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In the above sum rule we have con ned us to term s of lowest order in perturbation and 


operators of dm ension Jess than six. 


The correlator ofthe 1* current J%,, BJ) and its conjugate has the form 
2 
ai dxe“ "nnr ek Di= œ el): (29) 


We nd for i, (!) the denticalsum ukasthat for 9, (!), as t should be. 


The corresponding fom ula when the current Jo; 2 and Ji; 2 I Pd) and DI) are used 
instead of J%, 2 and J? p 2 are the follow ing 
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The correlators for the (1° ;2* ) doublet have the fom 


i ak) GQ +5=32(l); (32) 
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W hen the currents PÅ) and BJ) are used the sum mks for this doubkt are found to be 
Z 


P 1 lte pe 1 s ; 
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Here theh—G*item comes from thediagram (c) in Fig. 1 which vanishes forthe doubkt 
(om ees 


From (29), BJ) and BJ) we obtain curves shown in Fig. 2a,b,c. Forthe rst two cases 
stabilty window s exist when !, lies in the interval 2:4 to 2:9 GeV; forthe last one !.. les 
in the nterval 2:7 to 32 GeV . M posing usual criterum forthe upper and lower bounds 


of the Borel transform param eter T we found approxm ately the same stability w indow 
0:7{0:9GeV .The results for for three cass are the follow ing 


(=27)= 1:05 005 GW; (36) 


for the doublet (0* ;1* ) when the currents (20), a) w ithout the derivative is used. 


(1=2*)= 09 010 G&V; (37) 


for the sam e doublet when the currents [24), [29) with the derivative is used. 


B=2*)= 0:95 010 GW; (38) 


for the doublet (11 ;2*). 


The predicted m ass value for the doublet (0* ;1*) in (34) is som what di erent from 
that in B7). The value in Bq) is in agem ent with the result of [7] in which the same 
currents as DI) and (27) are used. Ba) predicts a mass value of the doublet (0* ;1*) 
slightly higher than them ass value ofthe doubkt (1* ;2* ) predicted in Ba) . On the other 
hand, B7 gives a mass value for the doublet (0* ;1") approxm ately equal or slightly 
lower than that for the doubkt (1* ;2*) in BJ) . Though (36) and (37) can be considered 
as consistent w thin the uncertainties, (37) indicates that the QCD sum me approach can 
incorporate amass value for (0° ;1* ) lower than that for (1* ;2* ). This is in agre% ent 
with quark m odel predictions in [11] and B]. By contrast, we have also derived the sum 
we forthe ground state doublet (0 ;1 ) with the cumentsw ith the extra derwative ið. 
W e found that there is no acceptable stability window for this sum we. 


Finally, we would like to address the follow ing points of this paper. The results from 


ourprevious work on B-S wave functions are used only for nding the appropriate currents. 


Once these currents are given, the calculations in QCD sum wis are independent of the 
consideration of the B-S equation. That the currents used are appropriate was checked 
by proving (11) with direct calculations in tem s of quark-glion elds nm HOET without 
referring to the B-S wave functions. A ssum ing duality and inserting intem ediate heavy 
meson states between J and JY in (11), this equation also strongly suggests equation (7) 
which is the direct statem ent that the currents create pure states w th de nite values of 
j,P and >}. 
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Figure C aptions 


Fig. 1. Feynman diagram s contributing to the sum we. There is no 
Fig (c) when the currents (20) (21) are used. 


Fig. 2. Dependence of on the Borelparam eter T for di erent values 
of the continuum threshold !.. (a) (0) Pr (0* ;1*) doublet with the 
interpolating currents Z) gå) and Rg RI used respectively. From 
top to bottom the curves correspond to (A) !. = 2:;2:6;24 GeV and 
(0) !.= 2:9;2:65;24 GeV . (C) Prdoubkt (1° ;2* ). From top to bottom 
the curves correspond to !. = 3:1;2:9;2:7 GeV respectively. 
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